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I.  INTRODUCTION 

This  report  describes  the  development  of  a unified  rocket  nozzle  computer 
program  that  Is  Intended  to  serve  as  an  aid  In  the  study  of  the  physical  and 
chemical  processes  that  control  rocket  exhaust  plume  dynamics.  Including  plume 
visibility. 

Low  smoke  propellants  produce  a small  mass  fraction  of  primary  smoke 
particles  at  the  nozzle  exit;  hence  condensation  In  the  plume  is  critical  to 
the  production  of  secondary  smoke  and  thus  plume  visibility.  Primary  smoke 
particles,  though,  play  an  important  role  in  plume  visibility  by  serving  as 
nucleatlon  sites  for  secondary  smoke,  and  therefore  a complete  determination 
of  plume  visibility  requires,  in  addition  to  equilibrium  condensation  boundary 
considerations,  an  accurate  description  of  nozzle  exit  plane  particle 
properties. 

To  improve  nozzle  exit  plane  predictive  capability,  a new  computer  code 
has  been  developed  by  incorporating  the  two-phase  transonic  flow  solution  of 
Kllegel  and  Nickerson*  (TSNIC)  into  the  AeroChem  FULLNOZ  code.“  This  develop- 
ment allows  the  use  of  chamber  initial  conditions  to  generate  a continuous 
nozzle  calculation,  thus  eliminating  the  assumptions  of  supersonic  start  line 
gas  and  particle  properties  previously  required  by  FULLNOZ. 

The  use  of  FULLNOZ  in  the  supersonic  regime  offers  powerful  calculational 
capabilities  by  simultaneously  treating  gas/partlcle  nonequilibrium,  non- 
equlllbrlum  chemistry,  diffusion  across  streamtubes,  and  (optional)  the  turbu- 
lent boundary  layer.  The  new  code  will  enable  the  accurate  determination  of 
particle  spatial  distributions,  velocity,  and  thermal  lags  at  the  nozzle  exit 
plane  and  provide  complete  input  information  for  subsequent  plume  calculations. 

The  following  sections  give  a brief  description  of  the  techniques  employed 
in  TSNIC  and  FULLNOZ.  A parametric  study  has  been  made  to  show  the  sensitivity 
of  the  solution  to  variations  in  (1)  supersonic  start  line  parameters  due  to 
transonic  approximations,  (ii)  particle  sizes  and  size  distributions  over  a range 


1.  Kliegel,  J.R.  and  Nickerson,  G.R. Axlsymmetrlc  Two-Phase  Perfect  Gas 
Performance  Program,"  TRW-02874-6006-R000,  April  1967. 

2.  Pergament,  H.S.  and  Thorpe,  R.D. , "A  Computer  Code  for  Fully-Coupled  Roc- 
ket Nozzle  Flows  (FULLNOZ),"  AeroChem  TP-322,  AFOSR-TR- 75-1563 , NTIS  AD 
A019  538,  April  1975. 
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applicable  to  low  signature  propellants,  and  (iii)  particle  drag  and  heat 
transfer  coefficients  over  the  applicable  range.  A complete  description  of 
program  input,  a sample  input,  and  a sample  output  are  given  in  Appendixes  A, 
B,  and  C,  respectively. 


II.  CODE  DESCRIPTION 


A.  TSNIC 

TSNIC  is  fully  described  in  Ref.  1.  A perturbation  method  is 
employed  to  attain  an  approximate  solution  to  the  subsonic  and  transonic  gov- 
erning equations  of  gas/particle  equilibrium  flow.  The  solution  includes 
velocity  and  thermal  lags  between  gaseous  and  inert  condensed  phases.  One  di- 
mensional sink  flow  is  assumed  in  the  converging  section  of  the  nozzle,  and 
gas /particle  properties  are  determined  on  the  sink  line.  Tlie  transonic  region 
is  divided  into  several  zones  and  perfect  gas  assumptions  are  used  to  obtain 
gas  properties.  Particle  trajectories  are  calculated  through  the  throat 
region  to  determine  particle  properties  on  the  supersonic  start  line.  The  gas/ 
particle  expansion,  including  nonequilibrium  effects  is  approximated  in  each 
zone  with  an  average  expansion  coefficient  (see  Fig.  1). 

TSNIC  uses  as  input  chamber  properties  (such  as  produced  by  standard 
equilibrium  codes),  particle  sizes,  distributions  and  mass  loading,  and  throat 
geometry  parameters.  Values  of  the  gasdynamic  state  and  particle  velocity  and 
thermal  lags  are  interpolated  on  the  start  line  for  a specified  number  of  points. 

B.  FULLNQZ 

The  FULLNOZ  code  is  fully  described  in  Ref.  2.  It  is  based  on  the 
MULTITUBE  method  of  Boynton’  which  uses  the  streamtube  method  to  Integrate  the 
hyperbolic  governing  equations  of  supersonic  flow.  In  this  method,  the  ellip- 
tic Navler-Stokes  equations  are  reduced  to  hyperbolic  form  by  assuming  that 
dlffuslonal  effects  along  streamlines  are  small  compared  to  diffusion  across 
streamlines.  This  assumption  is  very  good  for  rocket  nozzle  flows  and  enables 
one  to  solve  an  initial  value  problem  (where  a marching  procedure  can  be  used) 


3.  Boynton,  F.P.,  "The  MUI.TITUBE  Supersonic  Flow  Computer  Code,"  Genernl 
Dynamlcs/Convalr  GDC-DBB  67-003,  February  1967. 


4 


TP-368 


FIGURE  1 TSNIC/FULLNOZ  SCHEMATIC 


rather  than  the  more  difficult  boundary  value  problem.  The  gas  flow  equations, 
in  finite-difference  form,  are  solved  along  and  perpendicular  to  streamlines 
while  a full  continuum  particle  cloud  system  of  equations  is  incorporated  for 
the  condensed  phase. 

The  advantages  of  using  the  streamtube  method  over  the  method-of- 
characteristics  in  calculating  rocket  nozzle  flows  are:  (1)  species  diffusion, 
shear,  and  heat  transfer  normal  to  streamlines  are  easily  included;  (li)  chemi- 
cal reactions  or  internal  relaxations  are  easily  Incorporated  since  the 
calculation  follows  streamlines;  (iii)  bounding  surfaces  and  gradients  normal 
to  streamlines  are  treated  without  difficulty;  and  (Iv)  a wide  variety  of 
boundary  conditions  Including  mass  transfer,  shear,  and  heat  transfer  can  read- 
ily be  Incorporated. 

FULLNOZ  utilizes  the  following  information:  (1)  initial  gas  and 
particle  properties  in  the  supersonic  region  just  downstream  of  the  nozzle 
throat  (supplied  by  TSNIC) , (2)  the  nozzle  wall  contour,  (3)  a chemical  reac- 
tion mechanism  and  rate  coefficients,  (4)  physical  properties  of  the  particles, 
and  (5)  the  nozzle  wall  temperature  (for  the  boundary  layer  calculations). 
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With  the  above  input  data  the  marching  scheme  steps  from  one  orthogonal  sur- 
face to  the  next  (Fig.  1),  computing  gas  and  particle  properties  within  a 
specified  number  of  streamtubes.  The  (essentially)  continuous  particle  size 
distributions  are  . presented  by  up  to  a maximum  of  eight  discreet  particle 
sizes.  All  particles  cross  the  gas  streamlines,  but  the  lighter  particles 
follow  the  gas  streamlines  more  closely  than  the  heavy  particles.  As  depicted 
in  Fig.  1,  limiting  particle  streamlines  are  computed  for  each  particle  size. 

At  each  orthogonal  surface  (i.e. , each  integration  step)  the  code  calculates 
wall  shear  stress  and  heat  transfer,  boundary  layer  displacement  thickness, 
and  velocity  and  temperature  profiles.  The  shear  stress  and  heat  transfer  are 
coupled  to  the  main  nozzle  flow  via  their  effect  on  the  wall  streamtube  proper- 
ties. This  code  approximately  treats  the  effects  of  shocks  which  might 
originate  from  the  nozzle  wall  in  turning  the  flow  as  strong  compression  waves.* 
The  operation  of  the  combined  TSNIC/FULLNOZ  code  is  straightforward, 
but  does  require  a preliminary  equilibrium  calculation  to  generate  chamber  gas 
properties.  Because  TSNIC  does  not  treat  chemical  kinetics,  chemical  equilib- 
rium (or  interpolation  between  equilibrium  and  frozen  flow)  can  be  assumed  to 
extend  to  the  throat  region  to  obtain  the  composition  data  used  by  FULLNOZ. 

This  procedure  provides  acceptable  accuracy  in  many  situations  and  equilibrium 
codes  such  as  the  RPL-ISP  program*  are  especially  suitable  for  providing  such 
information.  If  greater  accuracy  is  required,  a separate  kinetic  calculation 
can  be  made  between  the  chamber  and  throat  using  a code  such  as  the  AeroChem 
nonequilibrium  streamline  code  (NEQSLINE) . ® A complete  description  of  program 
input  is  given  in  Appendix  A. 


* Nozzle  shocks  could  be  incorporated  into  FULLNOZ  in  a manner  similar  to 
that  employed  in  the  AIPP  Code*  to  detect  internal  plume  shocks,  but  this 
would  require  additional  programming. 

U.  Pergament,  H.S.  and  Kelly,  J.T.,  "A  Fully-Coupled  Underexpanded  Afterburn- 
ning  Rocket  Plume  Program  (The  AIPP  Code).  Part  II.  Program  User's  Manual," 
Final  Report,  AeroChem  TP-328,  AFRPL-TR-75-52 , NTIS  AD  A019  411,  December 
1975. 

5.  Selph,  C.  and  Hall,  R. , "Theoretical  ISP  Program,"  AFRPL,  Edwards,  CA,  1974. 

6.  Pergament,  H.S.  and  Mlkatarian,  R.R.,  "Prediction  of  Mlnutemnn  Exhaust 
Plume  Electrical  Properties,"  AeroChem  TP-281,  May  1973. 
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III.  CODE  TESTING 

In  order  to  determine  performance  characteristics  for  the  TSNIC/ 
FULLNOZ  code,  a parametric  set  of  calculations  was  made  to  test  for  the  sensi- 
tivity of  the  solution  to  variations  in:  (i)  start  line  properties,  including 
those  that  result  from  assumptions  concerning  throat  geometry  and  start  line 
location,  (ii)  particle  sizes  and  mass  loading  which  were  varied  over  a range 
expected  in  rocket  nozzle  flows,  and  (iii)  particle  drag  and  heat  transfer 
coefficients,  commensurate  with  the  expected  extremes  in  particle  shapes. 

The  parametric  calculations  were  made  for  a small  motor  using  the  N5 
double  base  propellant.  The  10. 5®  nozzle  has  a throat  area  of  0.72  cm“  and  an 
expansion  ratio  of  3.7.  To  save  computation  costs,  and  because  the  main  inter- 
est is  in  the  resultant  particle  properties  at  the  nozzle  exit,  the  comparisons 
were  made  with  no  kinetic  or  boundary  layer  calculations.  The  standard  case 
employed  the  following  values  of  key  variables:  (i)  mass  fraction  of  parti- 
cles, WPWGT  = 0.005;  (il)  ratio  of  axial  start  line  location  to  throat  radius, 
ZAX  = 0.7;  (iii)  particle  radii,  RP  = 0.8,  1.25  y;  (iv)  drag  and  heat  transfer 
coefficient  factors,  FFF  and  FFG  = 1.0;  (v)  nozzle  inlet  angle,  THID  = 10.5°; 
and  (vi)  ratio  of  throat  radius  of  curvature  to  throat  radius,  RRT  = 3.0. 

A.  START  LINE  PARAMETERS 

Key  input  specifications  to  the  TSNIC  code  that  affect  supersonic 
start  line  computations  are  the  throat  geometry  and  the  prescribed  location  of 
the  start  line.  For  many  cases,  the  details  of  the  throat  geometry  may  not  be 
known  and  estimated  values  must  be  input  to  the  program.  Hence,  calculations 
were  made  varying  the  throat  radius  of  curvature,  R^,,  and  the  converging  sec- 
tion inlet  angle  to  test  for  their  effect  on  start  line  and  ultimately,  exit 
plane  gas  and  particle  properties.  It  was  found  that  changing  the  inlet  angle 
from  10.5°  to  15°  produced  nearly  no  change  in  start  line  or  exit  properties, 
indicating  the  solution  is  insensitive  to  this  parameter.  There  is,  however, 
a restriction  in  the  selection  of  the  R^,  since  smaller  values  increase  the 
divergence  of  gas  and  particle  streamlines  at  the  throat,  complicating  the 
solution.  Reference  7 recommends  RRT  * Rc/R*  ^ 1*5  to  allow  the  solution  to 

7.  Coats,  D.E.,  et  al,  "A  Computer  Program  for  the  Prediction  of  Solid 
Propellant  Rocket  Motor  Performance,"  AFRPL-TR-75-36,  July  1975. 
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proceed  but  in  the  present  use  It  was  found  that  RRT  greater  than  2.0  was 

required.  It  is  believed  though,  that  this  restriction  does  not  alter  the 

adequacy  of  the  analysis  for  most  cases  of  interest. 

Another  restriction  in  the  use  of  TSNIC  is  in  the  prescription  of 

the  supersonic  sta..t  line  location.  The  axial  location  (downstream  of  the 

throat)  can  be  estimated  from  the  following  source  flow  approximation: 

0.63 

1 - cos(THIW) ' 


ZAX 


{ 


(RRT)sin(THIW) 


sin(THIW) 


') 


where  ZAX  =»  Z . /R*  is  the  start  line  location  and  THIW  is  the  nozzle  wall 
axis 

angle.  For  the  present  case  with  THIW  = 10.5®  and  RRT  = 3.0,  ZAX  = 0.7.  When 
ZAX  was  reduced  to  0.5,  TSNIC  failed;  when  a value  of  1.0  was  used  the  resul- 
tant start  line  was  geometrically  inconsistent  with  the  downstream  nozzle 
geometry  and  the  FULLNOZ  calculation  failed.  Other  considerations  in  the  use 
of  TSNIC  are  well  documented  in  Ref.  7. 


B.  PARTICLE  SIZES 

Early  studies®’®  of  primary  particle  sizes  in  rocket  motor  exhausts 
indicated  that  the  location  of  the  distribution  peaks  depends  on  factors  such 
as  motor  size  and  chamber  pressure.  For  small  motors  size  peaks  in  the  1 to 
3 u radius  range  were  expected.  More  recent  works'°>^‘  Indicate  that  the 
majority  of  particles  are  very  small  with  radii  less  than  0.1  u,  although  mass 
distribution  peaks  remain  in  the  1 to  3 y range. 

A series  of  calculations  was  planned  to  test  the  performance  of  the 
TSNIC/FULLNOZ  program  over  this  applicable  range  of  particle  sizes.  Unfortu- 
nately, though,  the  assumptions  employed  in  the  TSNIC  analysis  result  in 


8.  Smith,  P.W. , Delaney,  L. J. , and  Radke,  H.H. , "Summary  Results  of  Particle 
Size  Measurements,"  in  Proceedings  of  the  AFRPL  Two-Phase  Flow  Conference 
1967,  AFRPL-TR-67-223,  Vol.  I.,  p.  264. 

9.  Rochelle,  W.C.,  "Review  of  Thermal  Radiation  from  Liquid  and  Solid  Propel- 
lant Rocket  Exhausts,"  NASA  Marshall  Space  Flight  Center,  Huntsville,  AL, 
NASA-TM-X-53579,  February  1967. 

10.  Dawborn,  R.  and  Klnslow,  M. , "Studies  of  the  Exhaust  Products  from  Solid 
Propellant  Rocket  Motors,"  Arnold  Air  Force  Station,  Tullnhomn,  TN,  AEDC- 
TR-76-49,  September  1976. 

11.  Varsl,  G. , "Summary  of  Particulate  Measurements,"  NASA,  Atmospheric  Effects 
Working  Group  Meeting,  Vandenberg  AFB,  October  27-28,  1976,  p.  96. 
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serious  computational  difficulties  when  particles  with  radii  less  than  1 p are 
encountered.  Hence,  the  lower  limit  of  the  size  rnnp.e  tested  was  limited  to 
0.8  y.  The  upper  limit  of  2.5  y was  used  as  representative  of  the  small  motor. 

Comparisons  of  gas  and  particle  velocities  and  temperature  are  shown 
in  Figs.  2 and  3 for  particles  of  radii  0.8,  1.25  and  2.50  y.  A particle  load- 
ing of  0.5%  was  used.  Because  of  the  low  particle  loading,  the  particles  have 
a negligible  effect  on  the  gas  properties  and  the  profiles  are  uniform  at  the 
nozzle  exit.  The  particles  do  not,  however,  completely  follow  the  gas  stream- 
lines, as  evident  in  Fig.  4.  Tlie  angle  (j)  Is  the  flow  angle  of  the  particles 
relative  to  that  of  the  gas.  For  most  of  the  flow  regime,  the  particle  tra- 
jectories diverge  more  than  that  of  the  gas,  except  near  the  wall.  The  gas/ 
particle  divergence  Is,  in  all  areas,  quite  small. 

Increased  effects  of  the  gas/partlcle  Interaction  are  seen  when 
particle  loadings  are  Increased  to  30%  (characteristic  of  highly  aluminized 
propellants).  Velocity  and  temperature  profiles  at  the  nozzle  exit,  shown  in 
Figs.  5 and  6,  indicate  that  the  presence  of  particles  does  affect  the  exit 
plane  gas  properties.  In  these  figures,  case  A1  represents  a calculation  with 
particles  of  0.8  and  1.25  y radii,  and  case  A2,  1.25  and  2.50  y. 

Other  applications  of  the  FULLNOZ  code  to  variations  in  particle  size 
are  described  in  Ref.  12  for  a large  bell-shaped  nozzle. 

C.  PARTICLE  DRAG  AND  HEAT  TRANSFER 

The  momentum  and  energy  exchange  between  the  spherical  particle  and 
gas  flows  is  described  completely  in  Ref.  2.  In  order  to  include  nonspherical 
particles,  alterations  to  the  particle  and  heat  transfer  coefficients  are 
necessary.  Although  liquid  condensates  would  normally  be  spherical,  solidifi- 
cation and  agglomeration  of  certain  species  could  form  nonspherical  particles. 
Dynamic  shape  factors,  which  relate  the  drag  of  an  arbitrarily  shaped  particle 
to  that  of  a sphere,  are  available  in  Ref.  13.  For  the  extreme  shapes  of  pro- 
late and  oblate  ellipsoids,  the  increase  in  drag  is  within  50%  of  the  spherical 


12.  Thorpe,  R.D. , Pergament,  H.S.,  and  Hwang,  B.C.,  "NOx  Deposition  in  the 
Stratosphere  by  the  Space  Shuttle  Solid  Rocket  Motors,"  JANNAF  9th  Plume 
lechnology  Meeting.  CPIA  Publ.  277  (Applied  Physics  Lab.,  Johns  Hopkins 
fJnlv. , Silver  Spring,  April  1976)  pp.  317-352. 

13.  Fuchs,  N.A. , The  Mechanics  of  Aerosols  (Pergamon  Press,  New  York,  1964). 
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FIGURE  2 NOZZLE  EXIT  PLANE  GAS  AND  PARTICLE  VELOCITIES 
Particle  mass  fraction  = 0.5%. 
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FIGURE  3 NOZZLE  EXIT  PLANE  GAS  AND  PARTICLE  TEMPERATtrRES 
Particle  mass  fraction  ■ 0.5%. 
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FIGURE  4 DIFFERENCE  BETWEEN  PARTICLE  AND  GAS  TRAJECTORIES  AT 

NOZZLE  EXIT  PLANE 


drag.  Since  definitive  information  of  the  effect  of  particle  shape  on  heat 
transfer  was  not  available,  these  coefficients  were  assumed  to  be  proportional 
to  those  for  particle  drag.  The  results  of  this  calculation  are  shown  In  Figs. 
7 and  8 for  1.23  p particles.  An  increase  in  the  drag  to  the  upper  limit 
yields  a 10%  decrease  in  both  particle  velocity  and  temperature.  In  order 
to  more  fully  describe  the  sensitivity  of  the  analysis  and  to  account  for  un- 
certainties in  particle  data,  the  case  of  a decrease  in  these  coefficients  is 
also  shown. 
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FIGURE  5 THE  EFFECT  OF  PARTICLE  MASS  LOADING  ON  EXIT  PLANE  GAS 
AND  PARTICLE  VELOCITIES 

Mass  loading  Cor  standard  case  Is  0.5%;  for  cases  A1  and  A2,  30%. 


TEMPERATURE,  K VELOCITY,  cm/sac  (*IO®) 
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FIGURE  7 THE  EFFECT  OF  VARIATIONS  IN  DRAG  AND  HEAT  TRANSFER 
COEFFICIENTS  ON  PARTICLE  VELOCITY 

Particle  mass  loading  0.5%. 
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FIGURE  8 THE  EFFECT  OF  VARIATIONS  IN  DRAG  AND  HEAT  TRANSFER 
COEFFICIENTS  ON  PARTICLE  TEMPERATURE 

Particle  mass  loading  = 0.5%. 
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IV.  CONCLUSIONS 


The  Incorporation  of  the  Kllegel  and  Nickerson  transonic,  two-phase  flow 
code  into  the  AeroChem  nozzle  code  has  resulted  In  a powerful  tool  for  the 
calculation  of  nozzle  exit  plane  gas  and  particle  properties.  Although  some 
care  Is  required  to  properly  specify  throat  geometrical  parameters,  and  parti- 
cle radii  are  limited  to  ^1.0  p,  the  code  is  applicable  to  a wide  range  of 
conditions.  The  calculated  gas  and  particle  properties  at  the  nozzle  exit  for 
various  particle  sizes  and  mass  loadings  display  expected  behavior  for  conical 
nozzle  flows.  The  major  effect  on  particle  concentrations  at  the  exit  is  the 
size  distribution;  the  effects  of  uncertainties  in  the  particle/gas  interac- 
tion, resulting  from  variations  in  particle  shapes,  is  small. 

V.  RECOMMENDATIONS 

The  TSNIC/FULLNOZ  code  is  an  ideal  basic  tool  for  the  study  of  fundamental 
processes  that  control  the  distribution  of  primary  smoke  particle  concentra- 
tions at  the  nozzle  exit.  There  la  experimental  evidence  that  primary  smoke 
Is  strongly  affected  by  chamber  ballistic  properties  such  as  pressure,  burn 
rate,  mass  flow,  and  propellant  grain  size.“*  However,  there  are  little  or 
no  primary  smoke  results  described  in  terms  of  important  nozzle  processes  such 
as  high  temperature  condensation.  A condensation  model  applicable  to  nozzle 
flows  would  treat  the  smallest  particles  as  molecular  groups  which  would 
diffuse  with  the  gas.  These  groups  would  be  treated  as  particles  after  attain- 
ing a predetermined  size.  An  adequate  treatment  of  this  problem  requires  a 
transonic  analysis  with  the  capability  of  treating  chemical  kinetics  and  a 
wider  range  of  particle  sizes.  Particle  growth  mechanisms  such  as  coagula- 
tion might  Justifiably  be  Ignored  since  the  rates  of  these  processes  may  be 
too  slow  to  have  a marked  effect  in  the  nozzle. 

It  Is  recommended  then,  that  further  efforts  in  describing  primary  smoke 
be  directed  toward  the  modeling  of  high  temperature  condensation,  within  the 
framework  of  a unified  nozzle  model  such  as  the  TSNIC/FULLNOZ  code. 


14.  Placzek,  D.W. , "Effect  of  Interior  Ballistic  Parameters  on  Rocket  Exhaust 
Smoke,"  Rohm  and  Haas  TR-S-282,  November  1970. 
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APPENDIX  A 

PROGRAM  INPUT  AND  OPERATION 

The  TSNIC/FULLNOZ  program  was  coded  for  CDC  6600  or  7600  computers.  It 
is  run  In  two  job  steps  with  separate  Input  for  TSNIC  and  FULLNOZ.  Transfer 
of  start  line  properties  Is  made  from  TSNIC  to  FULLNOZ  via  scratch  flies:  the 
gas  properties  on  unit  2 and  particle  properties  on  unit  3. 

I . TSNIC  INPUT 

All  Input  to  TSNIC  Is  made  through  NAMELIST/DATA/ . The  Input  Items  are 


divided  Into  three  groups:  (1)  chamber  products  of  combustion, 

data,  and  (3)  Inlet  and  throat  parameters.  For  some  Input  Items 
assumed  by  the  program  and  usually  need  not  be  supplied. 

A.  PROPELLANT  DATA 

Item  Input  Quantity 

(2)  particle 

, values  are 

Units 

CAPN 

viscosity  temperature  exponent,  specific 
constant  pressure 

heat  at 

none 

GAMMA 

specific  heat  ratio,  C /C 

2p  8v 

none 

GMGO 

chamber  gas  viscosity  coefficient 

Ibs/ft  sec 

PC 

chamber  pressure 

psla 

PR 

Prandtl  number 

none 

TGO 

chamber  temperature 

'R 

WMOLG  molecular  weight  of  gas 

B.  PARTICLE  DATA 

none 

CPL 

particle  heat  capacity  (T  > T ) 

” "m 

• 

cal/mole  “K 

CPS 

particle  heat  capacity  (Tp  <:  T ) 

cal/mole  ®K 

D(J) 

the  diameter  of  each  of  n particles  Is  to 
so  that  D(l)  < D(2)...  < D(N)  set  D(N+I) 

17 
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Item 

Input  Quantity 

Units 

HPl. 

liquid  particle  enthalpy  (Tp  “ T ) 

kcal/mole 

HPS 

solid  particle  enthalpy  (Tp  - Tp  ) 

kcal/mole 

ITOT 

number  of  particle  groups,  N.  N ^ 8 Is  required 

none 

SMP 

particle  density 

Ibs/ft’ 

TPM 

particle  solidification  temperature 

“R 

WPWGT 

ratio  of  particle  to  gas  weight  flow 

none 

WPWT(J) 

particle  weight  flow  fractions  corresponding 
each  of  the  above  particle  diameters,  D(J) 

to 

none 

C.  INLET  AND  THROAT  PARAMETERS* 

Item 

Input  Quantity 

Units 

Assumed 

Value 

DZI 

particle  trajectory  integration  step  size 

none 

0.002 

DZMIN 

inlet  step  size  parameter 

none 

0.002 

NILP 

number  of  Initial  line  points 

none 

15 

RRT 

ratio  of  throat  radius  of  curvature  to 
throat  radius.  A value  R^.  2 is  required. 

none 

RT 

throat  radius 

ft 

THFD 

faring  angle  (no  faring  for  THFD  5»  THID) 

degrees 

5 

THID 

inlet  angle 

degrees 

THIW 

Intersection  of  Initial  line  and  wall 

degrees 

12. 

THJD 

angle  defining  the  zone  farthest 
downstream 

degrees 

9. 

ZAX 

ratio  of  Intersection  of  Initial  line 
and  axis  to  throat  radius 

none 

ZI 

number  of  upstream  zones 

none 

3 

ZJ 

number  of  downstream  zones 

none 

2 

* See  Fig. 

A-1 
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FIGURE  A-1  INLET  AND  THROAT  PARAMETERS 
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1 1 . FULLNOZ  INPUT 

Input  for  the  FULLNOZ  portion  of  the  program  is  given  below. 
Card  Fortran 


No. 

Cols. 

Name 

Description 

Format 

1 

1-80 

ID 

Run  identification 

20A4 

2 

1-10 

NSEC 

Maximum  run  time  (sec);  when  run  time 

1. 10 

reaches  NSEC  nozzle  properties  at  last 
orthogonal  surface  are  punched  and  can 
be  used  to  continue  the  calculation. 


(Not  operational  in  present  version  of 
code.) 

11-22  DXLSS  Approximate  distance  (cm)  along  axis  E12.5 

between  orthogonal  surface  print  sta- 
tions. For  control  of  print  increment 
by  number  of  integration  steps  (KP, 

Card  4,  Cols.  11-15)  instead  of  axial 
distance,  make  DXLSS  larger  than  XLMAX 
(below) . 

23-34  XLMAX  Maximum  distance  (cm)  along  axis  for  E12.5 

which  calculations  will  be  made.  XLMAX 
should  be  set  somewhat  larger  than  axial 
distance  to  nozzle  exit  plane  in  order 
to  complete  calculation  at  last  wall 
point . 

3 1-5  IBUGSH  Debug  printout  index;  J5 

0 - No  debug  printout 

1 - Extensive  printout  for  debugging 

6-10  ITURB  Turbulent  flow  Indicator;  15 

0 - Invlscid 

3 - Constant  values  of  turbulent  mixing 
parameters,  p,  Pr  and  Le  are  input 
on  Card  9. 

4 6-10  NN  Number  of  terms  for  Cp  polynomial  curve  15 

fit  plus  one;  for  curve  fits  supplied 
with  program,  NN  = 6.  If  fewer  coeffi- 
cients are  used  for  a given  species  NN 
is  left  at  6 and  zeroes  are  input  for 
the  missing  coefficients  (see  Card  7). 

11-15  KP  Output  control;  number  of  orthogonal 

surfaces  between  print  stations.  To 
control  print  via  axial  distance  (DXLSS, 

Card  2,  Cols.  11-22)  instead  of  orthogo- 
nal surfaces  set  KP  larger  than  LPLANE 
(Card  4,  Cols.  36-40). 
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Card 

No. 

Cols. 

Fortran 

Name 

Description 

Format 

4 

16-20 

MMAX 

Number  of  points  needed  to  describe  shape 
of  inner  boundary;  always  set  equal  to  2 
for  axis  boundary. 

15 

21-25 

NMAX 

Number  of  points  required  to  define  nozzle 
wall  contour  (50  maximum) 

15 

26-30 

NDS 

Total  number  of  gas  species  in  flow 
(25  maximum) 

15 

31-35 

NITER 

Maximum  number  of  iterations  allowed  for 
iterative  solutions  (e.g.,  calculation 
of  streamtube  properties  for  variable  y) 
Recommended  value:  NITER  • 50 

15 

36-40 

LPLANE 

Maximum  number  of  integration  steps  for 
entire  calculation 

15 

41-45 

IKINE 

Number  of  chemical  reactions  (40  maximum) 

15 

NOTE: 

IF  NO  BOUNDARY  LAYER  CALCULATIONS  ARE  TO  BE  MADE  CARD  5 
MAY  BE  LEFT  BLANK,  BUT  MUST  BE  INCLUDED 

5 

1-  5 

TWALL 

Nozzle  wall  temperature  (“K)  (assumed 
constant) 

F5.0 

8 

IBLFLG 

Boundary  layer  property  printout 

Indicator; 

0 - Print  Reynolds  no.,  thickness  and 

friction  parameters 

1 - Print  above  plus  velocity  and 

temperature  profiles 

11 

11 

IBL 

Boundary  layer  calculation  Indicator; 

0 - No  boundary  layer  calculation 

1 - Calculate  boundary  layer 

11 

6 

1-  9 

ALPHAH 

Factor  that  multiplies  maximum  stable 
step  size;  Recommended  value: 

ALPHAH  - 0.8 

F9.4 

10-18 

EPSLON 

Amount  by  which  streamtube  Mach  numbers 
must  exceed  one  in  order  for  the  calcu- 
lation to  continue.  Recommended  value: 
EPSLON  =0.01 

F9.4 

19-27 

TOL 

Convergence  tolerance  for  Iterative 

F9.4 

solutions;  Recommended  value: 
TOL  - 1 X 10"* 
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Card 

No. 

Cols. 

Fortran 

Name 

Description 

Format 

6 

28-36 

DELTA 

Metric  exponent*; 

0 r Two-dimensional  flow 

1 - Axially  symmetric  flow 

F9.4 

37-45 

ATOL 

Maximum  allowable  fractional  change  in 
streamtube  area  per  step.  Recommended 
value:  ATOL  = 0.1 

F9.4 

NOTE:  CARD  GROUP  7 DEFINES  THE  THERMODYNAMIC  DATA  FOR  EACH  SPECIES. 
THERE  ARE  2 CARDS  PER  SPECIES  AND  NDS  SPECIES 


7.1.1 

1-13 

14-26 

27-39 

40-52 

A(l)\ 

A(2)  1 
A(3)  f 
A(4)  J 

Li 

La  Specific  heat  polynomial  constants 

La  for  first  species,  (cal/mole  °K) 

L.. 

E13.5 

7.1.2 

1-13 

A(5) 

La,  Specific  heat  polynomial  constant 
for  first  species 

E13.5 

14-26 

A(6) 

Le,  Enthalpy  constant  of  integration 
for  first  species 

E13.5 

^ AH“  - 1^’*  Cp  dT^  (kcal/mole) 

• 

• 

27-39 

CS(1) 

La,  Entropy  constant  of  integration 
for  first  species  (cal/mole  °K) 

• 

7.NDS.1 
7. NDS. 2 

Repeat  thermodynamic  data  for  NDS  species 

NOTE: 

CARD  GROUP 
TRANSPORT 

8 IDENTIFIES  EACH  SPECIES  AND  SPECIFIES  VARIOUS 
PROPERTIES.  MUST  BE  IN  SAME  ORDER  AS  CARD  GROUP  7. 

8.1 

1-  4 

IDENT(l) 

Species  name;  first  species  (the  remain- 
ing data  on  Card  8.1  also  apply  to  the 
first  species) 

A4 

13-24 

MUO(l) 

Viscosity  at  reference  temperature^ 

E12.4 

(g/cm-sec);  only  used  when  flow  contains 
particles.  FOR  NO  PARTICLES,  MUO  MUST  = 0 
FOR  ALL  SPECIES 


* The  particle  conservation  equations  are  written  only  for  axially  symmetric 
flow.  Thus  two-dimensional  solutions  can  only  be  obtained  for  flows  with- 
out particles. 

^ Values  for  common  gases  can  be  found  in  most  physics  and  chemistry 
handbooks. 
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Card 

No. 

Cols. 

Fortran 

Name 

Description 

Format 

8.1 

25-36 

TO(l) 

Reference  temperature  for  viscosity  (“K) 

E12.4 

37-48 

OMEGA (1) 

Exponent  describing  viscosity/temperature 
relation;  y ® T*^ 

E12.4 

49-60 

PR(1) 

Reciprocal  of  species  Prandtl  number 

E12.4 

61-72 

SC(1) 

Reciprocal  of  species  Schmidt  number 

E12.4 

• 

• 

73-80 

MW(1) 

Species  molecular  weight 

E8.2 

• 

8.NDS 

Repeat  for  each  species 

NOTE: 

CARD  9 IS  NOT  NEEDED  IF  ITURB  * 0 (CARD  3 COLS.  6-10) 

9 

1-10 

TLE 

Turbulent  Lewis  number  (constant) 

E10.3 

11-20 

21-30 

TPR 

EDDYK 

Turbulent  Prandtl  number  (constant) 

Eddy  viscosity,  g/cm-sec  (constant) 

E10.3 

E10.3 

10.1 

1-10 

C(l) 

Mass  fraction  of  first  species 

E10.3 

11-20 

C(2) 

Mass  fraction  of  second  species 

• 

E10.3 

71-80 

C(8) 

• 

Mass  fraction  of  eighth  species 

E10.3 

10.2 

Continue  on  subsequent  cards;  up  to  25 
species 

8E10.3 

NOTE: 

CARD  GROUP 
WALL  POINT 

11  DEFINES  THE  NOZZLE  WALL  CONTOUR.  FIRST 
(AT  START  LINE)  IS  TRANSFERRED  FROM  TSNIC. 

11.1 

1-12 

XB(2) 

Axial  position  (cm)  of  second  point  along 
nozzle  wall 

E12.4 

• 

13-24 

RB(2) 

• 

Radial  position  (cm)  of  second  point  along 
nozzle  wall 

E12.4 

• 

• 

11(NMAX-1)1-12 

• 

• 

XB(NMAX) 

Final  axial  position  of  nozzle  wall 

E12.4 

RB(NMAX) 

Final  radial  position  of  nozzle  wall 

E12.4 

NOTE: 

CARDS  12  AND  13  ARE  INCLUDED  ONLY  IF  PARTICLES  ARE 

IN  FLOW. 

12 

1-10 

FFF 

Factor  to  adjust  partlcle/gas  drag 
coefficient 

E10.3 

11-20 

FFG 

Factor  to  adjust  particle/gas  heat 
transfer  coefficient 

E10.3 

21-30 

CL 

Liquid  particle  specific  heat  (cal/g-“K) 

E10.3 
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Card 

No. 

Cols. 

Fortran 

Name 

Description 

Format 

12 

31-40 

CS 

Solid  particle  specific  heat  (cal/g-'K) 

E10.3 

i^l-50 

HTRAN 

Particle  heat  of  solidification  (cal/g) 

E10.3 

51-60 

WT 

Particle  molecular  weight 

E10.3 

13 

1-10 

RHSS 

Particle  density  (g/cm®) 

E10.3 

31-40 

TPS 

Particle  solidification  temperature  (®K) 

E10.3 

NOTE; 

THE  FOLLOWING  CARDS  CONTAIN  THE  REACTION  MECHANISM  AND 
RATE  COEFFICIENTS.  USE  ONLY  IF  IKINE  (CARD  4,  COLS. 
41-45)  IS  GREATER  THAN  0. 

14.1 

1-  4 

IZD(l) 

Species  A 

A4 

7 

+ sign 

8-11 

IZD(2) 

Species  B 

A4 

14 

+ sign 

15-20 

Blank  or  M 

A4 

21 

= sign 

22-25 

IZD{3) 

Species  C 

A4 

28 

+ sign  (if  needed) 

29-32 

IZD(4) 

Species  D 

A4 

35 

+ sign  (if  needed) 

36-39 

IZD(5) 

Species  E 

A4 

49-50 

IRR 

Reaction  type* 

12 

51 

IRT 

Rate  coefficient*  type 

11 

52-59 

RC(1) 

Pre-exponential  factor,  A 
(cm-molecule-sec  units) 

E8.2 

60-63 

RC(2) 

Temperature  exponent,  N 

F4.1 

• 

• 

• 

6U-12 

RC(3) 

Activation  energy,  B (cal/mole) 

F9.1 

14.IKINE  Same  as  above  for  IKINE  reaction 


* See  Appendix  A,  Section  III  for  reaction  and  rate  coefficient  types.  See 
Ref.  2 for  a complete  explanation  of  the  reaction  mechanism. 
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HI.  REACT1QN_AND  RATE  COEFFICIENT  TYPES 


Ten  possible 

reaction  types  are  Included 

In 

the 

program: 

Reaction  Type 

(1) 

A 

B + M 2 

C 

+ 

D 

(2) 

A 

+ 

B + M S 

C 

+ 

M 

(3) 

A 

+ 

B Z 

C 

+ 

D 

+ 

E 

(4) 

A 

+ 

B Z 

C 

(5) 

A 

+ 

M z 

C 

+ 

I) 

M 

(6) 

A 

+ 

B 

C 

•f 

D 

(7) 

A 

+ 

B + M -► 

C 

+ 

M 

(8) 

A 

+ 

B -► 

C 

+ 

D 

+ 

E 

(9) 

A 

+ 

B •+ 

C 

(10) 

A 

+ 

M 

C 

•f 

D 

+ 

M 

The  forward 

rate 

coefficient,  k-. 

is 

input 

to 

the  code  as  one 

following  eight 

types 

Rate  Coefficient 

Type 

* 

(1) 

•^f 

- 

A 

(2) 

'^f 

AT"* 

(3) 

*‘f 

AT"* 

(4) 

‘'f 

AT**/* 

(5) 

A exp(B/RT) 

(6) 

‘‘f 

AT-*  exp(B/RT) 

(7) 

*^f 

AT-*/* 

(8) 

AT**  exp(B/RT) 

* Rate  coefficient  data  for  typical  rocket  nozzle  reaction*  may  be  found, 
e.g. , in  Ref.  IS. 

15.  Jensen,  D.E.  and  Jones,  G.A. , "Gas-Phase  Reaction  Rate  Coefficients  for 
Rocketry  Applications,"  Rocket  Propulsion  Establishment  Technical  Report 
No.  71/9,  October  1971. 
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The  equilibrium  constant,  Kp,  ig  determined  from 

In  Kp  - - AG/RT 

where  the  Gibbs  free  energy,  AG,  for  individual  reactions  is  computed  from  the 
input  thermodynamic  data. 
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APPENDIX  B 
SAMPLE  INPUT 

The  following  1h  Che  input  data  file  for  the  ntendard  case.  The  first 
file  la  Che  TSNIC  namelist  input,  and  the  second  file  is  the  FULLNOZ  input. 


•DATA  CPfA  CPS-A.  72.  OAf«tA-l  292.  HPL«3  I66.  HP8-3  027. 

pc-iaio.  • unoLO-aa.aa.  8nr«70D  . too-aiaa.  . tph-iocio  . itoi *2. 

WPyOT*  009.  UPWT«  9.  .9.  NT*.  0t97.  THJI)-7  . TNIU-IO  9.  THID-IO.  5. 


- W9  - 2. 
7 0 


79  IN.  NOTON 


90  fOOO 


TIN1C/FN02  SAm.C  TC6T  C 
1000  I.  40 

0 0 

A 9999  9 2 

0 

0 8 01  . 0001  1 0 0 1 
. A44009A2  £♦!  10979207  £♦!-. 440A909AC«0  3A200991 
94040940C’2-.  204439048  2.  94410327  8 2 
.907413318  1 139924048  2-0409440981  .219004978  1 

- 309470438  I-  944390728  2 943214998  2 

. 49403429  8 1-  233012488-3.  99447920  8-4-.  910729498-9 

- 104473238-4  90410499  8 2.  33404494  8 2 
90810249  8 1.  17399379  8 I-.  329342108 

. 79177941  8-1.  79212009  8 1.  90930994  8 
.40994491  8 1 1341219248  1-  149007928 
44414242  8-1-  170090298  1.  30443271  8 
.40311014  8 1.47049790  8 1-942912498 
40421192  8-1-  999939348  2.  91307249  8 
.42443449  8 1 19320092  8 1-.  441902748 
. 91092740  8-2-.  192424008  1.  92024774  8 


8-1 


0 22313947  8-1 
2 

0 00740497  8-2 
2 

0 49007994  8-1 
2 

0.  37003139  8-1 
2 


CO  1 
CO  2 
C02  1 
C02  2 

H I 
H 2 


OH 

Oii 

H2 

H2 

H20 

H20 

N2 

N2 


CO 

139  0 

8-4300. 

0 79 

1 

4 

26  0 

coo 

139  0 

8-4300. 

0.  79  . 

1 

4 

44  0 

H 

03.  9 

8-4300 

0 79 

1 4 

1 

4 

1 00 

OH 

129.  9 

8-4300. 

0 79 

4 

17  0 

H2 

03  9 

8-4300. 

0.  79 

1 

4 

2 00 

MOO 

129  9 

8-4300. 

0 79 

1 

4 

10  0 

lO 

140  0 

8-4300. 

0.  79 

4 

20  0 

9.908-1  1 708- 

4 498-4 

1 148-9 

1 908-2 

1 338-1 

1 318-1 

4 909  1 

7440 

1.  1 

. 0334 

0324 

9 907 

207  2 

11  341 

400 
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PMTICLC  KCVNOLM  NO 
PARTICLE  TEMPERATURE 


